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Abstract. We introduce and investigate a distance-type measure of non-Gaussianity 
based on the quantum fidelity. This new measure can readily be evaluated for all pure 
states and mixed states that are diagonal in the Fock basis. In particular, for an M- 
photon added thermal state, an analysis of the Bures degree of non-Gaussianity is made 
in comparison with two previous measures built with the Hilbert-Schmidt metric and 
the relative entropy. We obtain a compact analytic formula for the Hilbert-Schmidt 
non-Gaussianity measure and find a good consistency of the three examined measures. 
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1. Introduction 

Originally, non-Gaussian states were studied in quantum optics due to some non- 
classical properties: photon antibunching, quadrature or amplitude-squared squeezing, 
oscillations of the photon-number distribution. More generally, non-classicality was 
defined as the non-existence of the P representation as a well-behaved function. A 
review of the efforts made on these lines of research can be found in Ref.[T]. In the 
pure state case there exists a connection between non-classicality and the non-Gaussian 
charater of the density operator. Indeed, Cahill [2] proved that the only pure states that 
are classical are the coherent ones: all other classical states are mixtures. Therefore, all 
pure non-Gaussian states are non-classical. On the other hand, according to Hudson's 
theorem [3], all pure states with negative Wigner function are non-Gaussian. Signatures 
of non-classicality could be thus identified through the negativity of the Wigner function. 
It was shown that this holds for mixed non-Gaussian states as well [1]. Interest in the 
non-Gaussian states has recently emerged in quantum information processing. It was 
realized that non-Gaussian resources and operations could be more performant in some 
quantum protocols such as teleportation [HI El [7] and cloning [8]. To understand to what 
extent non-Gaussianity could be a resource in such cases, some distance-type measures 
of this property were proposed P [10] following one of the patterns 

5[p\r^ d^{p,PG) or 5[p\:= d{p,pG), (1.1) 

where pc is the Gaussian state having the same average displacement and covariance 
matrix as the given state p. In Refs. [U [TOl [11] , Genoni et al. employed as distance d 
the Hilbert-Schmidt metric and the relative entropy: they defined the Hilbert-Schmidt 
measure of non-Gaussianity [9], 

2 Tr(p2 

as well as the relative entropy of non-Gaussianity [TO] , 

8re[p\ ■■= S{p\pg) ■■= Tr(plnp) - Tr(plnpG)- (1-3) 

Interestingly, non-Gaussianity in terms of relative entropy [lO] was experimentally 
measured for single-photon added coherent states pi2]- Another approach to non- 
Gaussianity was based on Q function and lead to a measure expressed by the difference 
between the classical (Wehrl) entropies of the Gaussian state pc and the given non- 
Gaussian state p [13] . 

In this paper we introduce a measure of non-Gaussianity of a single-mode state p in 
terms of its Bures distance to the Gaussian state pc having the same first- and second- 
order moments of the canonical quadrature operators. In other words, our definition 
is of the type (11. ip and uses a well-known metric related to the fidelity between two 
quantum states . Fidelity-based metrics have proven to be fruitful in quantum optics 
and quantum information as measures of nonclassicality [I5] , entanglement [161 113 IB] , 



SHslp]:="-mS^^ (1-2) 
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and polarization [13 IMl Ell 1221 123] • We also intend to compare the three above- 
mentioned distance-type measures in analyzing the non-Gaussianity for a definite class 
of one-mode states. 

The plan of our paper is as follows. In Sec. 2 we introduce and examine the Bures 
degree of non-Gaussianity. We insist on its advantageous form for pure states and for 
mixed ones that are diagonal in the Fock basis. Section 3 investigates an interesting 
mixed state of this type which is important for experiments: a thermal state with M 
added photons. We here give a compact analytic form of 5hS) while 5re and dp are 
expressed in terms of two series which have to be summed numerically. Section 4 is 
devoted to a discussion of our numerical results, making a comparison between the 
three above-mentioned non-Gaussianity measures and analyzing their consistency. 



2. Bures measure of non-Gaussianity 

Following the definition (11. ip one can take advantage of the distinguishability properties 
possessed by the distance d in order to get reliable values for non-Gaussianity. For 
further convenience we write down the previously defined degrees [9l [in]. On the one 
hand, Eq. (II. 2p gives an easily computable expression: 



. Tr[(p-pc)^] _l 
^Hs[p\ - 2Tr(/)2) - 2 



Tr(p^) - 2Tr(/3Gp) 
Tr(p2) 



(2.1) 



On the other hand, despite its not being a true distance, the relative entropy is 
acceptable and used as a measure of distinguishability between two quantum states. 
Moreover, recall that the relative entropy of non-Gaussianity, Eq. (II. 3p . reduces to 

6re[p\=S{pg)-S{p), (2.2) 

where S{p) := — Tr(plnp) is the von Neumann entropy of the state p. Invariance 
properties of the degrees (12. ip and (12. 2p were discussed in detail in Refs.[9l UHl ITT] . 
We now define a fidelity-based degree of non-Gaussianity 

h[p] ■■= \ dlip, Pg) = 1 - VnP,PG)- (2.3) 

The explicit expression of the fidelity between the states pi and p2 was written by 
Uhlmann [IT] : 

^(pi,p2) = {Tr[(v^p2V^)V^]}'. (2.4) 

As seen in Eq. (12. 3p . fidelity is tightly related to the Bures metric ds introduced in 
Ref.[2l] on mathematical grounds. Several general properties of our definition (12. 3p are 
listed below as arising from well-known beneficial features of the fidelity |25] . 

(i) The property of the fidelity to vary between and 1 implies: 

Sf[p] = 0, iff p is Gaussian, (2.5) 

< 6f[p] < 1, otherwise. (2.6) 
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(ii) If at least one of the states is pure, Eq. f l2.4p reduces to the usual transition 
probability Tr(pi/32). Correspondingly, the Bures degree of non-Gaussianity (12. 3p 
of a pure state 1^1') (^^l is 

<5f[|^)(^|] = 1- v/(^Ipg|^). (2.7) 

(iii) As shown in Refs.[9l [TOl [11], when JJ are the unitary operators of the metaplectic 
representation on the Hilbert space of states, then p' = IJpW =^ {p')g = 
UpcU^ and, therefore, according to the invariance of the fidelity under unitary 
transformations we obtain the identity 

5p[UpU^] = 5f[p]. (2.8) 

It follows that Sf[p] does not depend on one- mode squeezing and displacement 
operations. 

(iv) The multiplicativity property of the fidelity has an interesting consequence on our 
definition (12. 3p for a two-mode product state pi ® p2- Indeed, if p2 is a Gaussian 
state, we get J^{pi ® P2, {pi)g ® h) = ^{Pi-, (Pi)g) and therefore 

^f[pi<^P2] = ^f[pi]- (2.9) 

(v) For commuting density operators, [pi, P2\ = 0, Eq. (12. 4 p simphfies to 

^(pi,P2) = [Tr(pfpf)]2. (2.10) 

Let us now remark that the properties (12. 5p and (12. 6p justify the interpretation of 5f[p\ 
as a degree of non-Gaussianity. At the same time, properties (12. 5p . (12. 8p . and (12. 9 p 
of 5f[p\ are shared by the non-Gaussianity measures (12. ip and (12. 2p as well [H]. Note 
that we do not discuss here the evolution of the non-Gaussianity of a state under a 
completely positive map which is expected to be a monotonic one in the cases of the 
relative-entropy- and fidelity-based degrees [25] . 

How well these measures discriminate between quantum states in order to be 
declared good measures of non-Gaussiannity is a complicated question which was already 
invoked when discussing distance-type measures of non-classicality [15] or entanglement 
[T6] . It is desirable that, for a specific family of non-Gaussian states, any distance- 
type degree has a monotonic behaviour with respect to the continuous parameters 
defining the set of states. In the case of one-mode states, we adopt as a reasonable 
criterion to verify the appropriateness of the non-Gaussianity measures (12. ip . (12. 2p . 
and (12. 3p . their monotonic behaviour with respect to the average photon number of 
the state (A^). Another property that one could expect for the three measures is 
their consistency, namely, their quality to induce the same ordering of non-Gaussianity 
when considering a specific set of states. It was already shown that relative entropy 
and Hilbert-Schmidt measures display different ordering for Schrodinger cat-like states 
[TT] . However, conclusions on such important aspects of distance- type degrees of non- 
Gaussianity cannot be drawn in general, but only for special sets of states. This happens 
because obtaining compact analytic results is a task that requires diagonalization of the 
density operator p followed by the exact summation of the corresponding power series. 
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Evaluation of 5hs[p\ seems to be easier than that of both 5re[p\ and 5f[p\- As a matter 
of fact, Uhlmann's expression (12 ■4p is not easy to calculate even on finite-dimensional 
Hilbert spaces. We refer the reader to a recent paper of two of us [26] where the state of 
the art in evaluating fidelity in the continuous-variable settings is presented. However, 
there are some important sets of states for which we can get explicit and relevant results. 
In the following we concentrate on such two computable cases. 

First, for pure states, Eq. ( 12. 7p shows that 5p is state-dependent. This is equally 
true for 5hs'- 

5Hsmm] = \ [l + Tr(p^) -2(vI/|pg|v&)] • (2.11) 

Indeed, in order to evaluate bp and bus-, we need to determine the reference Gaussian 
density operator pc and its expectation value in the pure state |\E')(\E'|. By contrast, 
the entropic non-Gaussianity measure fl2.2l) of any pure state is a unique function of 
a single variable, namely, the determinant of the 2x2 covariance matrix of the state, 
A := det(V): 

(2.12) 

It is worth mentioning, however, that the Fock states are the only pure states for which 
both the Bures and Hilbert-Schmidt degrees of non-Gaussianity (12. 7p and (12. lip depend 
solely on the parameter A. Let us consider a number state |M)(M|. The associated 
Gaussian state pc is a thermal one with the mean occupancy (iV) = M . Equations 
(12. 7p and (12.111) give, respectively, the formulae: 



and 



^F\\M)m\ = ^-\\j^^^^,. (2.13) 
M + 1 

fel|M)(M|l = ^-^j^^^, (2.14) 

with VA = M + i. Note that Eq. f imp was already derived in Ref.[9]. Owing to 
the invariance property (12. Sp . the expression (I2.13P coincides with the Bures degree of 
non-Gaussianity of squeezed or/and displaced number states. 
Second, for any mixed Fock-diagonal state, 

oo oo 

P = $^P/IOai with 5^P« = 1, (2.15) 

the Gaussian reference state pc is a thermal state with the same mean occupancy 
(iV) = Zp/. We denote o := {N)/{{N) + 1) and write its spectral expansion: 

oo ^ 

PG = J]sJO(^l With si = -^-—aK (2.16) 

1=0 \^^ / + 1 
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The corresponding Hilbert-Schmidt and entropic non-Gaussianity measures were written 
in Refs.[9| [To]. For further use, we cast the Hilbert-Schmidt measure (12. ip into a shghtly 
modified form: 



5 



HS\ 



1 + 
1 

2 



1 + 



1 



(2.17) 



Here we have used the purity of the thermal state pc arising from Eq. (I2.16p . while 
Qp{y) '■= ^iPiy^ is the generating function of the photon-number distribution of the 
given state p. The relative entropy of non-Gaussianity, Eq. (II. 3p . simplifies to: 



5re[p\ = '^{Pi In Pi - si Ins/) 



oo 

J2pi Inp/ + {{N) + 1) ln((iV) + 1) - (iV) ln((iV)). 



1=0 



(2.18) 



1=0 



In the last line we have used the von Neumann entropy of a thermal state. In this 
special case we notice the commutation relation [p, pc] = 6, which allows us the use of 
Eq. fl2:T0|) to get 



Sl- 



(2.19) 



Note that some important mixed non-Gaussian states have the structure (I2.15P : phase- 
averaged coherent states and various excitations on a thermal state of the type p ~ 



Pth{(i^)^of^ ■ Here a and o) are the amplitude operators of the field mode. 



3. An example: Photon-added thermal states 

In general, the states with added photons are non-classical and non-Gaussian. We 
choose here to analyze an M-photon-added thermal state [271 128] as an interesting 
example of a Fock-diagonal state whose non-classicality was recently investigated in 
ingenious experiments [291 EOl [31] . Its density operator is 

1 



P 



{M) 



Here M is the number of added photons and pth is a thermal state whose mean number 
of photons is denoted by n: 



Pth = {l-x)^x^\l){l\ with x:-- 



n 



1=0 



n + 1 



(3.2) 



Accordingly, the density operator p^^^\ Eq. (13. ip . has the following eigenvalues: 

I 



PI ■■= (p(^))« 



M 



l-x)^+V-*^, (/ = 0,1,2,3,...). (3.3) 
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The mean occupancy is simply 

(iV) = n{M + 1) + M, (3.4) 
such that the photon-number probabihties of the associated thermal state read: 

si-{PG [(M + l)(n + l)ri' (^-0,1,2,3,...). (3.5) 

The generating function 



oo oo / , \ 

1=0 l=M \ / 

has a compact form: 

Hence, the Hilbert-Schmidt scalar product of the states p^*^^ and p^^-* is 



Tr 



^(M)^{A/)j =yp^s, = -^-^ . (3.8) 

" J ((AT) +n+ 1)^+1 ^ ' 

We have still to evaluate the purity of the state p*^*^^: 

oo oo / , \ *^ 

i=0 /=M V / 

A change of the summation index in Eq. (13.91) leads us to a closed-form result 
proportional to a Gauss hypergeometric function, Eq. flA.ip : 

Tr [{p^^'^f] = (1 - xy^^'+^\Fi{M + 1, M + 1; 1; x^). (3.10) 

By applying the linear transformation (]A.2p . we eventually get the purity as a function 
of the ratio x, in terms of a Legendre polynomial (1A.3I) : 

T.[(p<«.)^].(i^)""V„(i±f!). (3.11) 

Note that the Legendre polynomial in Eq. (13. lip is strictly positive because its argument 
is at least equal to 1. Insertion of Eqs. (13. 8 p and (13.110 into Eq. (12.10 yields the compact 
formula 

Af + l T 



M 



2 VMim 



4(iV)+2 ((A^) +n + 1)^+1 

(3.12) 

with the mean occupancy (A^) given by Eq. (13.40 . The situation is different for both 
the entropic and Bures non-Gaussianity measures. Making use of the photon-number 
probabilities (13. 3p and (13. 5p in Eqs. (I2.18P and (I2.19p . we established noncompact 
formulae for the relative entropy and the Bures degree of non-Gaussianity. Each of 
their expressions includes a power series which has to be summed numerically. We 
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further computed numerically these two expressions as one-parameter functions of a 
single variable for several values of the parameter. 

Long ago, Agarwal and Tara [27] examined the non-classicality of the state (13.11) 
by writing its non-positive P representation and Mandel's Q-factor. Non-Gaussianity 
of this state was recently evaluated in Ref . [13] by employing the Wehrl entropy-measure 
and found to be equal to the non-Gaussianity of the number state |M)(M|, being thus 
independent of the thermal mean occupancy n. This is a consequence of an invariance 
property of the Wehrl entropy under a uniform phase-space scaling of the Q function of 
the state. 

4. Discussion and conclusions 

On physical grounds, we expect that a good measure of non-Gaussianity has a monotonic 
behaviour with respect to the mean photon number (iV) and, in turn, to the parameters 
entering its expression. It is quite clear that the non-Gaussianity measures f l3.12p . 
f l2.18p . and fl2.19p depend on the thermal mean occupancy n, unlike the Wehrl-entropy 
measure [13]. Our analytic formula, Eq. f l3.12p . led us to accurate values for the Hilbert- 
Schmidt degree of non-Gaussianity. In Fig. 1 we plot the three distance-type measures 
as functions of the parameter x for several values of the number M of added photons. 
It is interesting that the three measures of non-Gaussianity 5hs^ ^re, and Sp decrease 
monotonically with x. We did not find any extrema of these functions in contrast with 
Figs. 3 and 4 in Ref.[l3]. 




XXX 



Figure 1. Dependence of the distance-type non-Gaussianity on the thermal parameter 
X for A/-photon-added thermal states with M = 1, 3, 5, 10 (from bottom to top). 

The variation of non-Gaussianity with the number M of added photons is shown in 
Fig. 2 for several values of the thermal mean occupancy. Besides showing a monotonic 
dependence on both parameters n and M, Figs. 1 and 2 seem to display a consistent 
relation between the three non-Gaussianity measures involved. To better outline this 
aspect and inspired by Ref. [H], we plot in Fig. 3 their mutual dependences when the 
parameter x varies on its domain x G [0, 1] at the same values of the number M of 
added photons as in Fig.l. We can see that consistency is not present for all values of 
the parameters, especially in the dependence Sp — Shs- 

To conclude, in this paper we have introduced the Bures degree of non-Gaussianity 
built with Uhlmann's fidelity between the given state and its associated Gaussian 
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M M M 

Figure 2. Dependence of the distance-type non-Gaussianity on the number of added 
photons. All the plots start from origin. We have used n ~ 0.1,1,2,5 (from top to 
bottom) 




<5hs ^hs 5re 



Figure 3. Mutual dependence of the distance-type non-Gaussianity measures (|3.12|) . 
([238]) and (|2l^ for M = 1, 3, 5, 10 (from top to bottom). We have used n e [0, 1]. 



one. We have then investigated the behaviour of three distance-type measures of non- 
Gaussianity for an M-photon-added thermal state as functions of the variables M and n. 
We have found adequate monotonic dependences of 5hsi ^re, and 6f on both parameters 
M and n. This is displayed by Figs. 1 (as a function of the thermal mean occupancy) 
and 2 (as a function of the number of added photons). Although very different as 
geometric significance, the three measures seem to give consistent results by inducing 
the same ordering of non-Gaussianity. Figure 3 shows a very good consistency between 
6 RE and 6hs (left plot) and between 6f and Sre (right plot). We also notice that the 
plots corresponding to different numbers M of added photons are very close for all 
mutual dependences. 
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Appendix A. Some useful formulae involving Gauss hypergeometric 
functions 

A Gauss hypergeometric function is the sum of the corresponding hypergeometric series, 
,F,{a,b;c;z):=f2^^^^, (|z| < 1), (A.l) 
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where (a)„ := T{a + n) /T{a) is Pochhammer's symbol. This definition is extended by 
analytic continuation |32]. Recall the linear transformation formula 



The Legendre polynomial of degree M can be expressed in terms of a Gauss 
hypergeometric function: 
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